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Abstract We address the occurrence of narrow planetary rings and some of their 
structural properties, in particular when the rings are shepherded. We consider the 
problem as Hamiltonian scattering of a large number of non-interacting massless 
point particles in an effective potential. Using the existence of stable motion in 
scattering regions in this set up, we describe a mechanism in phase space for the 
occurrence of narrow rings and some consequences in their structure. We illustrate 
our approach with three examples. We find eccentric narrow rings displaying sharp 
edges, variable width and the appearance of distinct ring components (strands) 
which are spatially organized and entangled (braids). We discuss the relevance of 
our approach for narrow planetary rings. 

Keywords Narrow planetary rings • scattering systems • strands • braids 



1 Introduction 

The accidental discovery of the ri ngs of Uranus in 1977 by stellar-occultation 
measurements (lElliot et al.. 1 1 19771) led to a renewed interest in planetary ring 
systems. The main reason was that the Uranian rings turned to be extremely 
different from those of Saturn: They are narrow, opaque, s harp-edged, inclined 
and e ccentric (jE lliot and Nicholso n. II 19841 : lEsposito.1l2002l) . To quote some fig- 
ures (iMurrav and Dermott. .i l999i) . the widest ring of Uranus, the e ring, is 20-96 
km wide with a nominal semi-major axis at 51,149 km; in comparison, the main 
rings of Saturn are a few thousand kilometers wide. The eccentricity of the e ring 
is 0.0079; Saturn's F ring has an eccentricity 0.0026. The structural properties 
have r aised a number of questions, most of which remain unanswered dEsposito. I 
120021 : ISicardvT1l2005h . For instance, an eccentric inclined narrow ring like the e 
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ring is expected to spread in rather short time scales, /niax = 6 x 10 years, which 
is "considerably smaller than the age of the solar system" (lEsposito. Il2002h . This 
is due to interparticle collisions, drag and differential precession. This suggests an 
efficient confinement mecha nism that mainta ins these structural properties of the 
ring over longer time scales dEsposito. Il2002h . 

To explain these structural features, new models were introduced where the 
confinement was induced by nearby moons. A mong these models we mention i n 
particular the shepherding model introduced by Goldreich and Tremaine. I (Il979l) . 
where two moons around the ring were proposed to bound it. Soon afterwards, 
the Pioneer and Voyager missions detected the shepherd moons around Saturn's F 
ring, and later the shepherds of the outermost e ring of Uranus. The shepherding 
confinement involves angular momentum transfer between the shepherd moons 
and the ring particles, self-gr avity and viscous damping due to interparticle col- 
lisions dBorderies et al.. ' 1983). While the full scenario for shepherding has not 
been fully understood (Esposito, 2002; Sicajdy, 2005), the presence of dissipa- 
tion seems to be essential. In addition, the formulation assumes the ring boundary 
located at a lower-order resonance. 

While the discoveries of the shepherd moons around the F Saturn ring and 
the e ring of Uranus represented some confirmation for the shepherd theory, other 
questions still remain unanswered. For i nstance, most Uranian rings ha ve no as- 
sociated shepherd moons around them dMurray and Thomason. 1ll990h. nor the 
narrow eccentric rings of Sat urn, which among o thers would provide an explana- 
tion for their sharp-edges (.Murray and Dermott. il999 ). Thus, either the shepherds 
are the re but are too s mall to be detected, or "some physics is yet to be under- 
stood" ('Sicardy, "2005^ Other problems are related with the effects of collisions, 
drag or keplerian differential velocity, which define too short life-times for the 
narrow rings, and cast doubts on the confinement mechanism which maintains the 
eccentricity and sharp edges of these rin gs. In addition, Saturn's F' ring turned out 
to hav e a very rich dynamical structure (ISmith et al.. Ill98ll 1 19821 : iMurray et al.. I 
Il997h : besides the non-zero eccentricity, it displays multiple components entan- 
gled in a complicated way, known as strands and braids, showing further puzzling 
features like kinks and clumps. 

Numerical simulations have investigated a variety of physical interactions, 
like the g ravitational perturbations of shepherd moons on circular and eccen- 
tric orbits dShowalter and Bums. II 19821: loiu liatti Winter et al., 2000), effects due 
to the action of embedded moonlets (Lissauer and Peale. 1986), and ring inter- 
particle collision effects ( Hanninen, 1993; Lewis and Stewart, 2000). The cen- 
tral questions investigated have been the formation of structure (strands, braids, 
clumps) and their short-term stability. While these studies are extremely valuable 
and have led to interesting predic tions, e.g. the formation of channels and strea m- 
ers dCriuliatti Winter et al.. l2000l) which have been recently observed dMurray et al.. I 
I2OO5I1 there is no self-consistent approach for the confinement of narrow rings and 
their radial and azimuthal structure. The F Saturn ring remains as the most fasci- 
nating example of these possibilities. 

The present paper addresses this point, namely, a self-consistent scenario for 
the occurrence of narrow rings and the appearance of structure. The main idea is to 
consider phase space regions where scattering dominates the dynamics in systems 
with some intrinsic rotation. Scattering provides a simple scenario where confine- 
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ment and escape are clearly distingushed. Existence of stable periodic orbits or tori 
in such scattering regions implies effective (over long times) dynamically bounded 
motion where otherwise there would be escaping trajectories. The intrinsic rota- 
tion, which is associated to an external potential due to the shepherds or by the 
full many-body problem, defines a mechanism for creating the ring: the trapped 
region, if it exists, rotates around the central planet, thus the confined particles 
form a ring. While at first sight it may appear a strange idea to focus on scattering 
dynamics, close approaches to shepherds illustrate this situation. In plain words, 
what does not escape to infinity (due to scattering) builds a ring. The bifurca- 
tion scenario yielding the stable periodic orbits or tori, and other considerations 
supported by theorems of Hamiltonian systems with many degrees of freedom, 
allow us to understand the narrowness and the appearance of strands and braids. 
Within this approach, which at this stage is mainly qualitative, we succeed in ex- 
plaining the confinement and also some of the questions described above, such 
as the occurrence of narrow eccentric rings with sharp edges, apse alignment, the 
appearance of different ring components (strands) and their spatial entanglement 
(braids). The narrow rings we obtain for the examples treated are stable although 
they may display dynamical evolution. 

To illustrate our approach we consider three examples. The first, a planar bil- 
liard system (a disk rotating on a circular or Kepler ellipse) with an arbitrary large 
number of non-interacting massless particles, is very simple, yet the relevant steps 
can be followed to a good extent analytically. So it is an interesting system that al- 
lows a pedagogical insight on many properties of the scattering approach to narrow 
rings. The second example in addition includes a central attractive \/r interaction; 
it helps to establish the connection to the full gravitational problem, showing the 
genericity and robustness of our approach. The third example is the circular re- 
stricted three-body problem with an extremely small mass parameter. With these 
examples we show the occurrence of rings, the appearance of radial and azimuthal 
structures, and the genericity and robustness of the whole approach. 

The paper is organized as follows: in the next section we describe our approach 
in rather general terms, avoiding formalities. In Sec.[3]we illustrate our theory with 
the first example, a scattering billiard system. The main steps in the construction 
of a (narrow) ring are described in detail, discussing separately circular and ellip- 
tic motion of the billiard disk. For the elliptic case we present new results related 
to the appearance of structure within the rings, in particular, strands and braids. In 
Sec. m we present results for the rotating disk with the addition of a central attrac- 
tive 1 /r interaction. This example permits us to understand the connection to the 
gravitational problem, which we exemplify with the circular restricted three-body 
problem in Sec.|5] Finally, in Sec.|6]we present conclusions and some outlook. 



2 Hamiltonian scattering approach to rings of non-interacting particles 

We begin by considering the planar circular restrict ed three-body pro blem (RTBP) 
of point particles with a very small mass parameter dSzebehelv. Ill967i) ; more com- 
plicated and realistic cases will be considered below. The larger mass (primary) 
is thus associated to the central planet, the smaller one (secondary) to a confining 
moon, and the massless particle to a particle of the ring, and therefore it does not 
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influence the motion of the other two. The motion of the massive bodies is cir- 
cular about their common center of mass. In an inertial frame with the origin at 
the central planet, the Hamiltonian describing the motion of the massless particle 
takes the general form 

H = \{Pl + Vo(|x|) + yi(|x- x^((/))|). (1) 

Here, X = {X,Y) denotes the position of the massless particle, Px and Py are the 
canonically conjugate momenta, Xj/((/)) denotes the position of the secondary, and 
yo(|X|) and Vi (|X — X£/((/))|) are the interaction potentials due to the primary and 
secondary, respectively. The Hamiltonian is explicitly time-dependent through the 
dependence of Vx upon (/) = (Ot, with (/) the angular position of the secondary de- 
fined with respect to the X-axis and ft) its constant angular velocity. By canonically 
changing to a rotating frame (with constant angular velocity ft)), the new Hamil- 
tonian becomes time-independent and is thus a constant of motion (cf. Eq. Q 
below). This is the well-known Jacobi constant (Szebehelv, 1967). 

We observe that the Hamiltonian ([l]) has a quite generic form. It consists of an 
integrable part Hq = H — Vi{X,Y,(j)) and a periodic time-dependent pe rturbation 
Vi (X , Y,0). Following the standard terminology of scattering theory dNewton.l 
1 19661) . we shall denote by Hq the free part of the interaction, which is com- 
pletely integrable. In the case of the RTBP, the potential Vq is proportional to 
the mass of the central planet, while V\ , the disturbing potential, is proportional 
to the mass parameter fi. For the known shepherd moons we have /j ^ 10^^ — 
10->° (fUral'skaya, 2003). It is thus legitimate to consider the motion as small 
corrections about a Kepler solutions (Vb), i.e., in a series expansion. In this case, 
we say that Vq defines the geodesies of the free interaction Hq (the conic solutions 
for the Kepler problem), while Vi acts as a small perturbation. 

A large variety of periodic orbits (in the rotating frame) exist for the circu- 
lar RTBP with small jj. (,Henon. .il997.) . Co nsidering the measured ecce ntricities 
of the planetary narrow rings (< 0.01, cf. iMurrav and Dermott. Ill999h . we are 
thus interested in those which have a very small eccentricity, and therefore are 
close to the circular or elliptic orbits of the two-body Kepler problem. In addition, 
these orbits must occur in regions in phase space where scattering motion and thus 
the possibility to escape to infinity do minates the dynamics, i.e., no zero-velocity 
curves may exist dBenet et al.. Ill998h . That is, we focus on phase space regions 
with no "potential barriers" that confine the motion of the massless particle. It can 
be shown that such periodic orbits i ndeed exist, and some are linearly s table in cer- 
tain intervals of the Jacobi constant (iHenon. lll997l : lBenet et al..1l998h . Moreover, 
these stable periodic orbits appear through saddle-center bifurcations by changing 
the value of the Jacobi integral (e.g., by decreasing it): a pair of new periodic or- 
bits is created at certain value of the Jacobi constant, one stable and one unstable 
(see Fig. [3^ below). By slowly changing the value of the Jacobi integra l, the local 
horseshoe changes with respect to the parameters that characterize it (Be net et al.. I 
[1998). Eventually, further reducing the Jacobi integral induces a period-doubling 
bifurcation cascade on the stable periodic orbit, which finally turns it into an un- 
stable one (cf. Fig. [3J5); further reduction yields locally a hyperbolic horseshoe. 
The key property here is that, in the intervals of the Jacobi integral where there 
is one stable periodic orbit, the manifolds of the unstable partner define barriers 
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which confine dynamically the motion of the massless particle (Fig. [3^). Such 
dynamical confinement actually takes place all around the central planet due to 
the intrinsic rotation induced by the motion of the shepherd. Hence, there is a re- 
gion in phase space where the solutions remain close to the stable periodic orbit, 
not escaping, irrespective of whether the actual motion is periodic, quasi-periodic 
or even chaotic. This behavior is generic for autonomous two degrees-of-freedom 
scattering Hamiltonian systems of the type considered; systems with more degrees 
of freedom are discussed below. 

We turn now to the occurrence of rings dBenet and Seligmarini2000h . Consider 
an arbitrarily large number (ensemble) of massless particles whose interaction 
is given by ([T|, but which do not interact between themselves. While the latter 
assumption is very restrictive, it is a first step in the theory. It essentially allows us 
to separate the many-body problem into a collection of independent one-particle 
problems, which is the first term in the series expansion of the full many-body 
Hamiltonian in terms of the mass of the ring particles. The initial conditions of 
such non-interacting particles are chosen close to the phase space location of the 
stable periodic solution described above. Yet, they differ slightly, in particular 
in the value of the Jacobi integral. We emphasize that we are not interested in 
collision orbits with a shepherd but in periodic orbits very close to it. 

We restrict now to one interval of Jacobi integrals where the whole bifurcation 
scenario takes place (from the creation of the pair of periodic orbits, until both 
define a hyperbolic structure). For a fixed value of the Jacobi integral, all particles 
whose initial conditions are outside the phase space region enclosed by the man- 
ifolds of the unstable periodic orbit rapidly escape to infinity along a scattering 
trajectory. In contrast, all those within the manifolds are dynamically confined, 
moving around the central planet. The distinction among these two cases is sharp. 
Thus, letting the (non-interacting) system evolve from an initial time f = a large 
number of particles escape, while the rest move close to the reference stable solu- 
tion. These statements hold for all values of the Jacobi constant within the consid- 
ered interval. The ring is then obtained, in the X-Y space for a given time, from 
the position of all the ring particles of the ensemble that remain close to the refer- 
ence periodic orbit, i.e., that are dynamically confined (Figs. |4l|7] and [8ll. Note that 
this argument does not require a priori any resonance condition. 

Some important properties of the ring are the following. First, the ring dis- 
plays sharp edges. This follows from considering only scattering regions in phase 
space. Indeed, the distinction between dynamically trapped particles and those 
with unbounded motion is rather clear typically after a few periods of revolution 
of the secondary mass. Yet, for very small values of the mass parameter certain 
subtleties arise, as we shall discuss in sections |4] and |5] due to the overwhelming 
attraction of the central planet. Second, the rings are in general eccentric. This can 
be understood by observing that the motion of each particle of the ring is close to 
the originating stable periodic orbit, which in turn is close to a keplerian ellipse 
of small eccentricity in the case of the RTBP, as we have chosen it. In the rotating 
frame, the latter has a well-defined periapse and apoapse. In fact, the ring resem- 
bles the shape of the periodic orbit pictured in the rotating frame. Notice that this 
is a collective property reflected in the X-Y space, in the sense that it is defined 
by the ensemble of ring particles which are trapped. Put differently, apse align- 
ment follows from the proximity and resemblance of the trapped trajectories to 
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the organizing reference orbit (stable periodic orbit). Third, the narrowness of the 
ring can be understood from the fact that the region in phase space (for a specific 
Jacobi integral) corresponding to dynamically bounded motion is typically quite 
small (Benet et al., 1998), as well as the interval of values of the Jacobi integral 
where the reference periodic orbit is stable. This is usually so also because of 
the presence of other moonlets. We emphasize that, despite the qualitative nature 
of these results, these properties are obs erved in real narrow planetary r ings, and 
some of them remain as open questions dEsposito. ll2002l : [Sicardv. Il2005h . 

Above, we have used the fact that the Hamiltonian (dJ is autonomous with 
two degrees of freedom. This permits us to compute directly the stability prop- 
erties of t he periodic orbits from the trace and the determinant of the linearized 
dynamics dYakubovich and Starzhinskiinil975h . Moreover, in this case the man- 
ifolds of the unstable periodic orbit indeed define a region of bounded motion if 
the partner is stable. One consequence is that, for the interval of Jacobi constant 
of interest, the probability of b eing trapped is small but strictly non-zero; this fol- 
lows from the KAM structure dAmold et al.. Ill988h for a two degree-of-freedom 
system (Fig. [3^). Furthermore, the scenario is structurally stable, i.e., there are 
no qualitative changes under small generic perturbations: the organizing centers 
(periodic orbits) will only be shifted in phase space. Therefore, oblateness of the 
planet or other effects can also be included in this scenario. The whole approach, 
so far, is robust. 

The question now is whether this approach can be extended to systems with 
more degrees of freedom, for instance, when other shepherd moons are also taken 
into account, as in the case of Saturn's F ring or the e ring of Uranus. For more 
degre es of freedom new phenomena appear, such as Arnold diffusion (lAmold. I 
119641) . Furthermore, the larger dimensionality of phase spa ce allows fo r other 
possibilities as regards the stability of the periodic orbits (Sk okos. 1l200lh . These 
properties seem to be a real constraint for extending our approach to more degrees 
of freedom. Yet, for many degrees of freedom, the same line of reasoning can be 
followed with some minor changes: instead of presenting the theory on the stable 
periodic o rbits, we do it on the stable quasi-periodic tori. According to a theorem 
of Jorba and Villanueva. , (1997a; 1997b), around lower-dimensional normally el- 
liptic tori there is a region of effective stability. Moreover, if the corresponding 
unstable solution is of center x center x ■ • • x center x saddle type, recent results 
by Wiggins and collaborators assert that the scattering dynamics can be under- 
stood similarly to the invariant construction of the auto nomous two-deg ree-of- 
freedom system dWiggins et al.. 11200 ihlUzer et al.. ll2002l ; IWaalkens et al.. II2004I) . 
These results suggest that, at least locally and in the sense of effective stability, 
i.e. for long enough times but perhaps not infinitely long, the underlying phase 
space structure is essentially the same as for two degrees of freedom. In this case, 
the construction of the ring can be carried out, even though we cannot assure in- 
finitely long confinement. A simple case where these statements can be checked 
is the planar elliptic RTBP for small but non-zero eccentricity. In this case the 
explicit quasi-periodic time-dependence of the potential Vi cannot be removed, 
leading to a phase space of larger dimensionality. As we shall illustrate below, a 
quasi-periodic perturbation Vi has important consequences for the ring structure: 
the ring turns out to have many components (see Fig. [5^) which are entangled in 
a complicate way (cf. Figs.|5})-d). 
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We finish this section by emphasizing that, despite the fact that we have con- 
sidered the case of a shepherd moon for concreteness, the intrinsic rotation defined 
by the circular orbit above, or any other rotating quasi-periodic motion, may be 
associated with other moons or with the effect of the whole many-body interact- 
ing problem. If this intrinsic rotation allows for the appearance of trapped motion 
embedded in scattering regions, which are somewhat localized, then narrow rings 
will occur. This is therefore a possible explanation, not requiring undiscovered 
shepherd moons, for the non-shepherded narrow rings of Uranus for instance. 



3 Example 1: A scattering billiard system 

We illustrate our theory first with a planar billiard system. While this system is a 
simple toy model, we avoid the numerical complications of dealing with \/r po- 
tentials where the relevant masses differ significantly (many orders of magnitude). 
Yet, the connection to the gravitational case can be established rigorously (SeclU. 
This toy model, as we shall show, has the great advantage that it provides insight 
into the scattering properties on which our approach relies, since the fact that there 
is no global confining potential enhances the possibility of escape. 

Consider the planar motion of a massless point particle bouncing off one cir- 
cular hard disk of radius d, which moves on a Kepler orbit. The center of the 
disk X(i((/)) describes a circular or elliptic Kepler orbit with one focus at the ori- 
gin (Fig. [T]l. We denote its radial position by R{<\>), with (/) the angular position 
along the Kepler orbit measured from the pericenter. Then, 7?((/))=fl(l— e^)/(l + 
ecos(0)) with c( = 1 is the semi-major axis and e the eccentricity. The Hamilto- 
nian expressed in an inertial frame can thus be written as Eq. ([T]|, with Vq = 0, V\ 
zero for |X — Xrf((/))p > cfi and infinite otherwise. We shall focus on the case of 
small or zero eccentricity. For non-zero e, due to the explicit time dependence the 
system has two and a half degrees of freedom and no constant of motion. 

By performing a canonical transformation of (dJ to appropriate rotating-pulsating 
coordinates ( Merlo. ..20041 the disk is left at rest at = (^,0), where ^ = a is 
the major semi-axis of the orbit of the center of the disk. The new Hamiltonian / 
is then given by 



For e = (circular orbit), J is time independent and thus a constant of motion (the 
Jacobi integral). We shall refer here to the numerical value of / in (O as the Jacobi 
integral, even for non-zero eccentricity where it is not a constant of motion. 

The dynamics of the billiard is straightforward. The particle moves freely on a 
rectilinear trajectory of constant velocity (the potential Vq = 0) until it encounters 
the disk; if it never collides then it escapes to infinity. When a collision takes place, 
the particle is specularly reflected with respect to the local (moving) frame of the 
disk. This defines the outgoing conditions at the disk, and the motion is rectilinear 
and uniform again. The precise result of a collision depends on the position on 
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Fig. 1 The scattering billiard: ^ is the position of the center of the disk, a denotes the position 
of the collision point on the disk, v is the magnitude of the outgoing velocity and Q defines its 
direction. The center of the disk moves on a Kepler ellipse, whose foci are shown on the X-axis 
as open circles (o), and its center by a cross (+). 



the disk where it occurs, the relative velocities, and for non-vanishing e, on the 
position along the ellipse ((/)). Collisions taking place on the front of the disk 
increase the (outgoing) kinetic energy of the particle, while collisions on the back 
reduce it. 

In the context of rings, we are interested in the particles which are dynami- 
cally trapped, the ring particles. In the present example, this can only be obtained 
through consecutive collisions with the disk. It is important to emphasize here that 
collisions with the disk do not model or imply physical collisions with the actual 
shepherds; they model interactions with the external potential that allows for con- 
finement. A convenient description of this situation is given by introducing the 
following quantities, defined at the collision point (Fig. [l]). The angle (/) defines 
the position of the center of the disk with respect to the X axis (inertial frame), 
and V is the magnitude of the velocity after the collision. The angle a denotes the 
angle formed by Xrf((/)) and the position vector of the collision point referred to 
the center of the disk. Finally, the angle Q is the outgoing direction of the velocity. 



3.1 Occurrence of rings: The circular case 

For e = 0, Eq. (|2]l is an autonomous two-degree-of-freedom Hamiltonian. The 
Hamiltonian flow defines a map, e.g., at the collision point with the disk. This 
map is open, in the sense that certain initial conditions may not have an associ- 
ated image, as due to escape. Trapped trajectories are associated with consecutive 
collisions with the disk. In particular, simple periodic orbits can be worked out 
analytically. The Jacobi integral is given by / = — v (7?sin + li sin(0 — a)) 
(ft) = 1). Due to the circular symmetry, all radial collisions (a = tt) conserve the 
kinetic energy. It is easy to calculate the condition that these orbits must fulfill to 
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Fig. 2 (a) Radial (consecutive collision) periodic orbits of the circular rotating billiard for 
n = 0, 1,2,3, in terms of J /{R — d)^. (b) Stability of the radial periodic orbits: the orbits are 
stable ii\TtDi^j \ < 2. 



be fixed points of the map (iMeyer etal.. Ill995h . 

Jn _ 2cos2 0+zi(/)sin(20) _ 2cos2 0(1 +4(/)tan0) 



(3) 



Here, n = 0, 1 , . . . denotes the number of full turns completed by the disk between 
consecutive collisions, and A<^ = {2n — 1)71 + 20 is the corresponding change in 
the angle 0. In Fig. [2^ we plot /„/ {R — d)^ for some values of n. However, the 
periodic orbits are a set of measure zero of the initial conditions. Therefore, we 
must also consider the behavior of nearby trajectories. To this end, we note that 
the characteristic curves in Fig. |2^ display one maximum and one minimum for 
each value of n. The occurrence of maxima and minima on these curves points 
to the appearance of consecutive collision periodic orbits through saddle-center 
bifurcations. 

The stability properties of a time-independent two-degree-of-freedom system 
can be obtained from the trace and the determinant of the linearized Poincare 
map at the fixed points ( Yakubovich and Starzhinskii, 1975). The transformation 
{ajc+i,Pk+i) = ^j{cCk,Pk) with pi^ = ~d ~ Rcosttk — v&\n[ak — Qk) defines an 
iterative canonical map. Here, G [0,27r] and pi_ G [— PmaxiPmax], where pmax ~ 
{2J+R^ + d^+2Rdcosay^^. Then, the information on the stability is completely 
contained in the trace o f the linearized matrix D3^j. One finds ( fMerlo. 1120041 : 
iBenet and Merlo. II2004I) 

^ . (40)2(1 -tan^e) -4(1 +Zi0tan0) 
TrD^y = 2+ ^ ^' ^ d/R 

Changes in the stability of the periodic orbits, i.e. bifurcations, occur at TrD^j — 
±2. Equatin g the r.h.s. of Eg. Q to 2 tu rns out to be equivalent to the condition 
d/„/d0 = dBenet and Seligman. Il2000 h, i.e., the condition defining the position 
of the minima and maxima of /„. The corresponding result for —2 yields a con- 
dition related to period-doubling bifurcations. Figure |2b shows the behavior of 
TrD^y in terms of 6. Note that a simple semi-analytical estimate of the width of 
the rings is obtained by projecting the curves TrD^ = ±2 onto the X-Y plane. 

Figure|2j) shows connected intervals in 6, and therefore in /, where the radial 
collision periodic orbits are strictly stable, i.e. \TrD^j\ < 2. For / within these 
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Fig. 3 Phase space structure of the scattering billiard on a circular orbit near the maximum of 
the n = curve (R=l and d = I /2). (a) The horseshoe displays an elliptic and a hyperbolic fixed 
point (// (R — d)^ = 0.29325). (b) Both fixed points are unstable and the dynamics is dominated 
by scattering events, although the horseshoe is not yet completely hyperbolic {J/{R — d)^ = 
0.29218). 
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Fig. 4 (a) Stable ring of non-interacting particles of the scattering billiard system on a circular 
orbit corresponding to the n = stable periodic orbit, (b) Zoom of a region of the ring. The black 
lines are the analytical estimates given by TrD^j = ±2. 



intervals, the phase space structure of the scattering system displays one elliptic 
fixed point, surrounded by typical KAM tori and some chaotic layers (Fig. [3^). 
This region in phase space is strictly bounded by the stable and unstable invariant 
manifolds of the companion hyperbolic fixed point. Particles whose initial con- 
ditions lie within this region display infinite consecutive collisions and are thus 
dynamically trapped. Initial conditions that do not belong to any region of trapped 
motion eventually escape to infinity. Consider the case n = for concreteness. 
After the appearance of the pair of periodic orbits (when TtD3^j = 2), by further 
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reducing the value of J, there occurs a value where IxDS^j = —2. There, the el- 
liptic fixed point becomes inverse hyperbolic, and a period-doubling bifurcation 
sets in: the region of trapped motion vanishes rapidly after further reducing J. In 
Fig. I3J5 we illustrate this case, plo tting the phase space stru cture corresponding 
to an incomplete Smale horseshoe (iRiickerl and Jung7lll994) . The dependence of 
TiD3^j = —2 upon d/R implies that the actual parameters of the billiard influence 
the width of the ring. Physically, this means that the parameters related with the 
motion and the mass ratios of the shepherd moons do indeed influence the width 
of the ring. Further reducing t he value of / yields a hyperb olic Smale horseshoe. 

We turn now to the rings (iBenet and Seligman. Il2000h . Consider an arbitrary 
large number of non-interacting particles. Their initial conditions are completely 
arbitrary except for the value of their Jacobi integral, which satisfies the condition 
\Tr D3^j \ < 2 for the n = curve. Then, the particles may collide with the disk at 
any value (j) along the circular orbit of the disk. Letting the system evolve, many 
particles escape after a few collisions; the remaining ones are dynamically trapped. 
Note that the measure of the latter set is strictly positive, that is, there is a non- 
vanishing probability of finding such initial conditions. The motion of the particles 
belonging to this set, the ring particles, may be periodic, quasi-periodic or chaotic. 
The ring is the pattern formed by these trapped particles in the X-Y plane (inertial 
frame). Figure |4] shows the pattern obtained after a few thousand revolutions of 
the disk, resulting from a large number of particles whose initial conditions are 
close to the maximum of the n = curve. The plot shows the {X,Y) position of 
all particles, in an inertial frame, which have not escaped at the time when the 
"photograph" was taken. 

In Fig.|4]we also plot the curves corresponding to TtD3^j = ±2, which serve 
as estimates of the ring shape. They do not coincide perfectly with the boundaries 
of the ring, since they are related to the appearance and destruction of the central 
elliptic fixed point and carry no information on the size of the island around the 
elliptic points. Yet, the nature of these curves and the scattering dynamics implies 
a sharp-edged ring. The narrowness of the ring is a consequence of the relatively 
small region, in phase space, of trapped motion; this also takes into account the 
(small) variations induced by changing /. Furthermore, the ring can be character- 
ized as eccentric, in the sense that there are two points, the periapse and apoapse, 
forming a line that intersects the origin, the line of apsides. Finally, we mention 
that the same ring may be obtained by an ensemble of particles whose initial con- 
ditions, while different, are very localized in a region in the X-Y plane. The time 
evolution in this case not only selects the particles of the ring, but also spreads 
them over the whole region of trapped motion (in phase space) on a short time 
scale. 

3.2 Rings in the case of elliptic motion 

In the case of elliptic motion of the disk, Eq. (O is explicitly time-dependent. The 
organizing centers in phase space are now stable tori, which are indeed found at 
least for small enough values of the eccentricity. Close to these stable tori, the 
structure of phase space is s omewhat sim i lar to that discussed above for circular 
motion; details are given in iBenet at al.. I (|2005). Hence, rings occur just as be- 
fore. For very small values of the eccentricity, the rings have essentially the same 
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Fig. 5 (a) Detail of the ring when the disk moves on an eccentric Kepler orbit (e = 0.00167). 
Whole ring represented in terms of r' (see text) for (b) t = 7^/10 (top right), (c) f = 37^/10 
(bottom left) and (d) t = T^I2 (bottom right). Here, = 2n is the period of the disk's orbit. 

Structure as in the circular case. Increasing e slightly rapidly shrinks the width of 
the ring. Eventually, we reach a value of the eccentricity for which the ring dis- 
plays gaps, i.e., regions of extremely low density of particles. This is illustrated 
in Fig.|5t. In this case (e = 0.00167), the gaps divide the ring into three compo- 
nents, which we call strands. They are a direct consequence of the stable tori that 
yield the ring, which arise from the quasi-periodic character of Vi in Eq. ([T]l. We 
observe that particles belonging to a given strand stay in the same strand, i.e., they 
do not move along other components of the ring. 

In Fig. |5j) we show a different representation of the whole ring for a specific 
time (given as a fraction of the disk's period T^). We plot the azimuthal position 
<i> of the ring particles, measured anti-clockwise from the collision point with the 
disk as origin, and r' = r — r^t, where r is their radial position and rjit an appropriate 
fit to the mean value. This figure shows clearly the intricate entanglement of the 
strands, reminiscent of the beautiful braids of Saturn's F ring. We emphasize that 
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the different strands and braids we observe are consequence of the fact that a stable 
quasi-periodic torus is now the organizing center of the trapped region in phase 
space. Figures |5};-d display the ring at other fractions of the period of the disk. 
Together these plots are a representation of the orbital motion of the ring, i.e., the 
dynamical behavior of the braids. 

We emphasize that the strands and braids we have obtained for this toy model, 
which display a dynamical evolution, are not short-lived structural properties. 
Thus, stable tori allow for the occurrence of structured narrow rings. Moreover, 
and as we shall show in the next sections, the ideas described here can be extended 
to more realistic situations. This is so because of the robustness of relying on con- 
siderations in phase space. The above description of the appearance of strands and 
braids should therefore hold in more realistic situations. 

4 Example 2: Rotating billiard with a central gravitational interaction 

In the previous section we illustrated, using a billiard system, our approach to the 
occurrence of narrow rings and how some observed structural properties arise. 
This example, while showing the robustness of the ideas of the approach, con- 
stitutes the simplest scattering system that is able to produce narrow rings, and 
as such has a pedagogical value. However, it lacks of the gravitational potential 
with the central planet. In this section we include a central attractive \/r interac- 
tion. As we shall show below, the present case establishes the connection to the 
gravitational restricted three-body problem, which we consider in Sec.|5] 

Consider the billiard system on a circular orbit with the addition of a cen- 
tral attractive gravitational interaction, i.e., Vq = — \/r between the (ring) particle 
and the central planet. In distinction to the pure billiard, where the free motion 
is rectilinear and uniform (the potential is Vo = 0), the particles now move along 
solutions of the Kepler two-body problem. Besides such (collisionless) solutions, 
there are also trajectories which collide with the disk. We focus on the latter, as 
they model rings close to shepherds; the former would trivially form a two-body 
keplerian ring, which is circular and wide. The collision solutions are constructed 
by segments of keplerian trajectories, glued together at the collision points on the 
disk, where the precise outcome of the collisions is taken into account. 

We proceed as before and compute the consecutive-collision periodic orbits 
in the rotating frame. To this end, we observe that the radial collisions conserve, 
as before, energy and angular momentum, and occur only at the innermost (r_ = 
R — d) or outermost (r+ = R + d) points of the disk. As mentioned above, in the 
present case we must also account for the ksplsnan free flight between the colli- 
sions. The present problem is thus equivalent to finding periodic orbits in a central 
gravitational field which collide consecutively with a moving point singularity on 
a circular orbit of radius r_ and r+, respectively. This singularity influences the 
motion of the third particle only through collisions. This is so as long as the tra- 
jectory avoids other collisions with the disk; if not the periodic orbit is simply 
destroyed. We note that this is ess entially the same as the circular RTBP with 
^ = 0, as defined bv lHenon. I ( Il968l) . 

These statements can be formulated quantitatively as follows. The solutions of 
the keplerian free flight are precisely conic sections. We consider for concreteness 
the elliptic motion. Using symmetry arguments, it is easy to show that consecutive 
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Fig. 6 Periodic orbits of the disk on a circular orbit with a central gravitational interaction for 
(a) d = 10^^ and, (b) d = 0. The insets show differences between the two cases. 

periodic orbits can be obtained from initial conditions (of the moving singularity 
and the third particle, see Henon, ,,1968) on the X-axis of the inertial frame. This 
axis is parallel to the major axis of the ellipse. The elliptic motion of the third 
particle, in an inertial frame, is parameterized by 

Xp = aoae{a2cosE - e^), (5) 

Yp = Ooaia2ae{l - ) sin£, (6) 

f = Ue {E - a2eeSinE). (7) 

Here, and are respectively the major semi-axis and eccentricity of the elliptic 
orbit, and E is the eccentric anomaly. The parameter Gq defines if the periapse 
lies on the positive (1) or negative (—1) Z-axes, ai determines if the trajectory is 
direct (1) or retrograde (—1) and whether the initial condition is at periapse (1) 
or apoapse (—1). Then, in order to obtain the periodic orbit one needs to match 
{Xp,Yp) with (r-tCOsf,r±sinf) and solve for and £«. With a = O0G2, O-q — 
sgn(sinT]) and p = ar;(l — a cos t cos t])'/^, the timing condition reads 

/ N r 9 / M Tsin"^77 
Fo(T,T];a,ar,) =p[T]p^-smT](cos77 -o-cost)] Jn-=^- (8) 

r± 

iHenon. I d 19681) computed first the solutions of ([8| for r+ = r_ {d = 0). 

In Fig.|6^ we show the consecutive collision periodic orbits for a disk of radius 
d = 10~^; in Fig. |6j) the case d = is plotted for comparison. At first sight, no 
difference can be observed. The ca se d = corre sponds precisely to the circular 
RTBP with mass parameter IJ. —0 dHenon. II 19681) . This naive remark has impor- 
tant consequences. We emphasize that we are interested not in the billiard problem 
with a central gravitational potential, which is only another exarnple, but o n the 
full gr avitational problem, with its simpler version, the RTBP. IHitzl and Henon. I 
(Il977h proved that the limit, in the case of vanishing /i, of the second species 
periodic orbits of the circular RTBP with small but non-zero /i are the critical 
points of the consecutive collision periodic orbits for jU = (Fig. 16})). This state- 
ment proves the validity of our approach for the gravitational circular RTBP, and 
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Fig. 7 (a) Ring generated by the Ci2(2) periodic orbit for the circular rotating billiard system 
with a 1 //• central interaction. The multiple components of the ring are a consequence of the 
structure of the organizing periodic orbit, which is closed after a few full turns of the disk, 
(b) Detail of a region of the ring. 



establishes the connection of our naive examples to the gravitational case. Using 
the arguments of Section |2l this validity can be extended to problems with other 
interactions as well as more degrees of freedom. 

Coming back to the billiard with the central gravitational interaction, we are 
interested in the regions of phase space where scattering may take place. In this 
example, scattering is possible for positive keplerian energy. In the 9 — J plane, 
scattering motion occurs if the Jacobi integral satisfies, in some 9 interval, the 
inequality / > max(/_,7+), with /+ = — (2ro)'^^sin 9 and 7- = — 1/ro (ro = r±, 
depending on the value of 9). We notice that for 7 > scattering is possible. 

To construct the rings, we proceed as before, requiring that the periodic orbits 
where scattering is achievable are, in addition, stable. Knowing the periodic or- 
bits it is easy to compute their stability. The same procedure then yields the rings, 
though some subtleties arise. These are related to the fact that the particle's energy 
must be positive for the particle to escape. The only way of changing its energy is 
through collisions with the disk (which is why non-colliding orbits form trivially 
a wide circular ring). However, the change in the energy may not be so large as 
to lead to escape after a single collision, so that repeated collisions are typically 
required for escape. Since between collisions the particle moves along a Kepler 
ellipse, and we are interested in a tiny disk, the time between consecutive colli- 
sions can be rather long. In a Poincare surface of section, defined when the particle 
crosses the x synodic axis (e.g, x vs. px), such long motions along a Kepler ellipse 
will be represented as something resembling a torus. The reason to consider such 
a Poincare surface of section is inspired by the RTBP However, these ave. fake 
tori, not true tori, since collisions will abruptly change the particle into a different 
one (if the energy is still negative) or lead to escape. Therefore, long integrations 
(about 10^ periods of the disk or even longer) are required in order to distinguish 
whether the trajectory is truly trapped, and not to "misguide the eye" with the fake 
tori described above. 
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In Fig. |7^ we pres ent the ring gen erated by the stable collision periodic orbit, 
denoted as Cni2) by iHenon. I (Il968h . for a disk of radius d = 0.0006. The disk 
radius was chosen as large as possible but such that the characteristic curve in 
Fig. displays a local maximum and therefore the saddle-center scenario that 
yields one stable and one unstable fixed points takes place. As observed in Fig.|7^, 
the ring displays multiple components through the inner loop, reminiscent of a 
strand structure, as shown in Fig.[5]for the rotating billiard on an elliptic orbit. The 
observed structure is inherited directly from the shape of the organizing periodic 
orbit in the synodic (rotating) frame; there, the organizing periodic orbit is closed 
after a few full turns of the disk. Figure |7j) shows an enlargement of a region of 
the ring. As expected, the rings are narrow and eccentric, and display sharp edges. 

5 Example 3: The circular restricted three-body problem 

We now present results of our approach to the occurrence of narrow planetary 
rings in the circular restricted three-body problem. We concentrate on the case 
where the mass of the secondary (shepherd) is 10^^, in units of the mass of the 
central planet, i.e. jx w 10~^. While this value is much larger than the actual masses 
of the known shepherd satellites, it yields qualitatively the same results with an 
enormous saving in computing time. 

The advantage of considering the rotating billiard disk with a central l/r po- 
tential is that we know the precise location of the consecutive collision periodic 
orbits. Moreover, its finite size and the curvature of the disk enhance the insta- 
bilities which finally lead to escape. In the case of the RTBP, even for very small 
mass parameters, this knowledge has to be obtained numerically, despite the fact 
that Henon's work for the case /i = gives a starting point. The numerics are nev- 
ertheless non-trivial, since the orbits of interest are quite close to collisions with 
the secondary, which not only complicates the numerics but requires a rather high 
level of accuracy. As we shall see below, we also require rather long integrations. 

We shall focus on non-collision orbits generated by the continuation of the 
Ci2(2) consecutive collision periodic orbit of Henon for jU = 0. This illustrates 
that indeed the rotating billiard with a l/r central attractive interaction provides 
the link between the naive rotating billiard sy s tem a nd the RTBP. Moreover, this 
also illustrates the result of i Hitzl and HenonTI d 19771) that the generating periodic 
orbits for jU = are the second species periodic orbits for finite but small values of 
/i . We emphasize that, the organizing centers and the trapped orbits we consider, 
and therefore the ring, do not display collisions with the secondary. We first use 
Fig. |6j) to locate approximately the Jacobi integral interval where the continua- 
tion of the Ci2(2) periodic orbit is stable. From here, numerically, we obtain the 
interval where the periodic orbit has a stable branch. We find that such a stable 
orbit exists for 0.400 <J< 0.453. Note that these values of the J acobi integral 
are much larger than the associated with the zero-velocity curves dBenet et al.. I 
[l998). Therefore, scattering dominates the dynamics unless the proposed mech- 
anism (existence of stable periodic orbits) creates regions in phase space of dy- 
namically trapped motion. That is, escaping trajectories largely fill phase space 
for these values of J, except when trapping occurs near a stable periodic orbit. Be- 
low, we shall focus in the dynamics close to, and including, one specific region of 
bounded motion. For the RTBP with very small values of /J the dynamics consists 
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Fig. 8 Ring generated by the Ci2(2) periodic orbit for the circular RTBP with fi 10^^. (a) Full 
view; (b) details of a specific region close to the shepherd. Note the dependence of the width 
with respect to the proximity with the shepherd. The shepherd and the central planet are denoted 
by dots. 



of almost-elliptic keplerian motion, that takes place during long times, followe d 
by an eventual close approach that finally leads to escape dBenet et al.. Ill998h . 
Therefore, the numerical simulations must distinguish the fake tori described in 
the previous section, which again requires long and accurate integrations. In our 
simulations we have checked for about 10^ periods of the secondary that no abrupt 
tori exchange takes place. 

Once the stability region (bounded motion) around the periodic orbit is located 
and also the approximate interval where it is stable, the construction of the ring is 
straightforward, following the ideas we have described above. Figures |8] show the 
whole ring and an enlargement of a region. All qualitative features described in 
the previous sections are observed once again: sharp edges and a narrow eccentric 
ring with many components. Again, these many components are inherited from 
the shape of the organizing periodic orbit in the synodic frame, which closes after 
a few full turns of the disk. Another important observation is that there is a vari- 
able width of the ring, with respect to the position closest to the sh epherd. This 
property is observed in real narrow planetary rings dEsposito. II2002I) . The expla- 
nation within our framework is simple: the region of trapped motion is very small 
close to the shepherd, for any value of 7; away from it, the two-body Kepler inter- 
action allows for longer excursions due to the different values of / involved. The 
projection onto the X-Y plane is wider as we move away from the neighborhood 
of the shepherd. 

We note that, for the interval of the Jacobi constant considered, other regions of 
dynamically trapped motion exist in phase space. In particular, there is a compar- 
atively large stable region associated with a period-one fixed point, and another 
one related to a period-two fixed point, both belonging to phase space regions 
where scattering dominates the dynamics. From our construction, these regions 
would also yield rings. Yet, these stable regions are extremely robust, in the sense 
that they are present for an extremely wide interval of J. Focusing on the period- 
one stable region, it essentially sweeps out every value of x in the synodic frame. 
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Therefore, except for a tiny gap which is opened by the shepherd, it yields a wide 
circular ring; this ring is of no interest in the context of the present work. It is the 
resuh of the dominance of the central l/r keplerian interaction. This result implies 
that the RTBP is not a realistic model, at least to describe the Uranus rings or the 
narrow rings of Saturn. We believe that a realistic model should not only consider 
the planet, a shepherd and the ring particle as we have done, but at least another 
body corresponding to one of the major moons of the planet. For Saturn's F ring 
as well as for the e ring of Uranus, perhaps a fifth body is also required to fully 
account for the two shepherd moons. 

6 Conclusions and outlook 

In this paper we have presented a general and self-consistent approach for the oc- 
currence of shepherded narrow planetary rings, and some consequences for their 
structural properties. Our approach is based on the local structure of phase space 
around stable (periodic or quasi-periodic) solutions in regions where scattering 
dominates the dynamics. The corresponding structure in phase space allows for a 
set of positive measure to exhibit dynamical trapped motion, i.e. a non-zero prob- 
ability to find such trapped motion; particles with initial conditions outside this 
small region escape to infinity along scattering trajectories. The intrinsic rotation 
provides a mechanism for dynamical confinement around the central planet. Here 
we have considered that such rotation is induced by the perturbations of a shep- 
herd moon, though our approach is more general since the intrinsic rotation can 
also be associated with interactions with other, perhaps more massive, moons, or 
with the combined evolution of the interacting many-body problem. A possible 
consequence is the explanation of the non- shepherded narrow rings of Uranus: the 
many-body problem is such that it defines several regions of trapped motion, one 
for each specific narrow ring. An ensemble of (non-interacting) particles will form 
a well-defined ring. Focusing in regions where scattering dominates the dynamics 
allows to sharply distinguish between trapped motion related to the occurrence of 
the ring, and open motion manifested as regions with essentially no material due 
to the escape of the particles. 

Within our approach, the rings will be sharp-edged, as a consequence of the 
scattering dynamics. The size of the stability region in the extended phase space 
determines how narrow the rings are. For comparatively small regions of trapped 
motion, the narrow ring will typically be eccentric, since the particles will be 
close to the organizing center, which in general displays such eccentricity. Semi- 
analytical estimates of the edges of the ring can be obtained by analyzing the 
bifurcation conditions that define the region of bounded motion in the extended 
phase space (within the sea of scattering trajectories) around the organizing cen- 
ter, at least for models with two-degrees of freedom. Typically, the bifurcation 
scenario that yields the regions of dynamically trapped motion corresponds to the 
saddle-center bifurcation, and by consequence of its fast development, the rings 
are narrow. We have illustrated our considerations with three distinct examples, 
aiming to show explicitly the robustness of the approach and its extension to and 
consequences for a case with many degrees of freedom, and have proved that our 
considerations hold also in the gravitational case as illustrated by the RTBP with- 
out physical collisions. More realistic interactions, e.g. including the oblateness 
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of the planetary, can be easily incorporated within our theory. While this does not 
yet explicitly prove the connection to real narrow planetary rings, which are com- 
plicated collisional many-body systems with gravitational and other interactions, 
such an extension seems to us feasible, since our ideas are based on general phase 
space considerations, and on results and theorems which apply to Hamiltonian 
systems with more than two degrees of freedom. 

We emphasize that our construction of the rings is self-consistent, in the sense 
that we do not introduce any "cell" or periodic boundary conditions on a seg- 
ment of the ring to make the problem amenable in any sense. Therefore, radial 
and azimuthal structures do appear. Moreover, the theoretical framework devel- 
oped is not restricted to short time scales, at least as long as the model consid- 
ered is believed to be applicable for longer time scales. It is in this scope that 
our results may be relevant and contribute to the understanding of real planetary 
narrow rings: the mechanism of confinement is very efficient. We have explic- 
itly showed which phase space structures allow for the occurrence of such nar- 
row rings. Moreover, we have showed how rings develop further radial and az- 
imuthal structure: rings may display variable width and many components which 
are entangled in a way that, at least qualitatively, resembles the observed braids 
of Saturn's F rings. Within our approach, the appearance of different strands and 
braids is not related to the ring interparticle co llisions, as proposed dHanninen. I 
Il993h . nor to the existence of other moonlets dLissauer and Peale.1ll986l) : they 
are related to the quasi-periodic character of the perturbing force of the rotat- 
ing interaction, as induced by the shepherd's explicitly eccentric motion. This 
non-circular motion makes the Hamiltonian system effectively have many de- 
grees of freedom, which complicates its analysis, but allows for other bifurca- 
tions as well as resonances. As far as we are aware, this is a new explanation 
for the occurrence of strands and braids, although it is somewhat in agreement 
with th e idea of associating periodi c clumps with the eccentricity of shepherd 
moons dShowalter and Bums! The strands and braids that we obtain dis- 

play a periodic dynamical evolution, where the relevant time scale is related to 
the time for repeated interactions, or close approaches, with the shepherd. In our 
opinion, this points out the importance of the quasi-periodic perturbation, i.e. the 
non-circular motio n of the in trinsic potential. This observation is in agreement 
with the results of Giuliatti Winter et al., (2000), although these authors focus 
on certain perturbations induced by the repeated interactions of a group of parti- 
cles with Prometheus within a cell of the ring. Whenever the organizing periodic 
orbit closes after more than one full turn of the secondary, the ring displays non- 
independent multiple components, reflecting the shape of the organizing periodic 
or quasi-periodic orbit in the rotating frame. Despite the fact that we have consid- 
ered only planar motion, and that this is again only a qualitative comparison, the 
generalization of these "higher periodic" orbits to a non-plana r case serves as an 
interpretation of re cent observational findings dCharnoz et al..' 1 ,20051 Moreover, 



preliminary results dBenet and Merlolin preparation!) show that our theory can also 



yield other observed properties of narrow rings, such as arcs qualitatively similar 
to those of Neptune (Namouni and_Pprco, 2002). 

Our ultimate goal is to provide quantitative comparisons with observed plane- 
tary rings. The results presented here are only qualitative due to the fact that, so far, 
the examples are analogous to restricted three-body like problems. Actually, the 
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most realistic example we considered is the circular (gravitational) RTBP. As we 
have discussed, non-collision phase-space structures of the RTBP with extremely 
small values of the mass parameter, which originate from the overwhelming dom- 
inance of the central planetary attraction, make such simple modeling definitely 
not realistic, since they yield two-body keplerian wide circular rings. We have 
proposed to consider the interaction of other planetary satellites and, perhaps too, 
other shepherds as a more realistic model. The shepherds can be included in a 
first step as zero-mass moonlets which break certain solutions of the many-body 
problem (planet inc luding oblateness, other m oons and a particle of the ring). Fol- 
lowing the ideas of iHitzl and Henon. I d 19771) for the RTBP with ^ = 0, periodic 
or quasi-periodic collision solutions with respect to the shepherds can be generat- 
ing orbits for the case of small masses of the shepherds. Considering the stability 
regions in the extended phase space when the shepherds have finite mass allows 
a straightforward construction of the associated rings. Work illustrating a many- 
body problem is in progress. 

Our considerations may also be used to extend the applicability of current 
models to longer time scales, once it has been showed that the ring particles con- 
sidered do belong to some regions of dynamically bounded motion. Conversely, 
they may allow to qualify certain rings or some of their structural properties as 
temporary phenomena, if a realistic model does not have dynamically trapped 
particles. An interesting question is whether chaotic motion of the shepherds may 
still yield long-lived rings within our framework, that is, whether they allow for 
bounded motion, and what the consequences are for their structure. This is related 
to recen t investigations where the rnotion o f Prometheus and Pand ora seems to be 
chaotic (iGoldreich and Rappaportni2003allbl : lRenner et al.. Il2005h . 

The strongest assumption within our approach is the non-interacting character 
of the ring particles. This essentially allows us to separate the many-body prob- 
lem into a collection of independent one-particle problems. This assumption must 
be taken as a first approximation in the series expansion of the full many-body 
Hamiltonian in terms of the mass of the ring particles. Here, the particle-particle 
interactions among the ring particles will be orders of magnitude smaller than 
other interactions, except maybe when they collide. It is precisely the inelastic 
collisions among the ring particles what we neglect. A possibility to include these 
collisions in a simple way is with a random-walk model, taking place in all the 
coordinates of the full phase space of the ring particles. Such a model could also 
be used to describe dissipative effects like drag, depending on the particular im- 
plementation. This random-walk approach would certainly introduce an erosion 
mechanism of the ring, i.e., the possibility of escape from the ring and a typical 
life-time; yet it will preserve the properties that we discuss, such as the sharp- 
edges, the eccentricity and the many components of the ring. This is because, 
even when collisions are included, once a particle is outside the region of trapped 
motion it escapes along a scattering trajectory. Work in this direction is also in 
progress. 
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